All-optical control of a solid-state spin using coherent dark states 
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The study of individual quantum systems in solids, for use as quantum bits (qubits) and 
probes of decoherence, requires protocols for their initialization, unitary manipulation, and read- 
out. In many solid-state quantum systems, these operations rely on disparate techniques that can 
vary widely depending on the particular qubit structure. One such qubit, the nitrogen-vacancy 
(NV) center spin in diamond, can be initialized and read out through its special spin-selective 
intersystem crossing, while microwave electron spin resonance (ESR) techniques provide unitary 
spin rotations. Instead, we demonstrate an alternative, fully optical approach to these control pro- 
tocols in an NV center that does not rely on its intersystem crossing. By tuning an NV center to 
an excited-state spin anticrossing at cryogenic temperatures, we use coherent population trapping 
and stimulated Raman techniques to realize initialization, readout, and unitary manipulation of 
a single spin. Each of these techniques can be performed directly along any arbitrarily-chosen 
quantum basis, removing the need for extra control steps to map the spin to and from a preferred 
basis. Combining these protocols, we perform measurements of the NV center's spin coherence, 
a demonstration of this full optical control. Consisting solely of optical pulses, these techniques 
enable control within a smaller footprint and within photonic networks. Likewise, this approach 
obviates the need for both ESR manipulation and spin addressability through the intersystem 
crossing. This method could therefore be applied to a wide range of potential solid-state qubits, 
including those which currently lack a means to be addressed. 

1 I 'o explore control of individual quantum states, our experiments exploit coherent dark reso- 



J- nances that occur in a basic quantum mechanical level configuration known as a lambda (A) 
system. This configuration, consisting of two lower energy states coherently coupled to a single 
excited state, has been observed in a wide array of systems including atoms HJ, trapped ions, di- 
amond nitrogen-vacancy (NV) centers [2, 3], quantum dots [4], superconducting phase qubits 0, 
and optomechanical resonators |6|. In trapped ions, A systems can additionally be exploited to 
drive stimulated Raman transitions providing unitary rotations of the qubit state (Zl HJ. This ver- 
satile structure also forms the framework for a variety of other important advances in quantum 
science such as electromagnetically induced transparency f9j, slow light flQ|, atomic clocks [11], 
laser cooling [12], and spin-photon entanglement IIT31 . 

Here, we use time-resolved methods and quantum state tomography to explore the dynamics of 
various optically driven processes within a solid-state A system (Fig. [l]^)- This allows us to demon- 




strate three all-optical quantum control [[Sj [TH [151 protocols for a single NV center: initialization, 
unitary rotation, and readout of its spin state. Our A system consists of two ground state spin 
sublevels coupled to a spin-composite excited state sublevel formed by tuning the excited states to 
an avoided level anticrossing. Driving transitions between the levels of our A system resonantly 
with appropriate coherent light fields (Fig. [TjA.) causes any initial mixed state to be purified! 16], 
or trapped, into a well-defined but selectable quantum superposition. This superposition is called 
the "dark state" since destructive interference from the driving fields causes the system not to be 
optically excited. This dissipative effect, known as coherent population trapping (CPT), allows us 
to initialize the precessing spin anywhere on the rotating-frame Bloch sphere, the geometric surface 
corresponding to all possible superposition states of the spin. Opposite the dark state on the Bloch 
sphere is a corresponding "bright state" which couples strongly to the optical fields. Together, these 
dark and bright states define a unique basis whose orientation within the rotating frame (SI Ap- 
pendix) is a function of the relative phase and amplitude of the two driving optical fields (Fig. [l]A). 
A complementary process allows us to read out the spin state within this selected basis because 
the resultant photoluminescence (PL) during the transient period of the CPT interaction is pro- 
portional to the spin's projection along the bright state. Furthermore, detuning the driving fields 
from resonance within the A system produces unitary rotations of the spin state about a chosen 
dark/bright-state axis, a dispersive technique that is a product of stimulated Raman transitions 
(SRT). Thus, this A system approach allows spin initialization, readout, and rotation schemes to all 
function within a fully mutable basis. 

A A System in the NV Center 

The negatively charged NV center consists of a substitutional nitrogen atom adjacent to a lattice 
vacancy within a diamond crystal flT7|l . Its millisecond-scale coherence times!18J are exceptional for 
a solid-state system, and coherent optical transitions enable important applications in quantum op- 
tics and quantum information processing such as spin-light coherence! 19 1 and entanglement! 13 1, 
single-shot readout[20], coupling to photonic cavities |2Tl |22 | , two-photon interference! 23. 24J, im- 
plementations of quantum games!25j, and photon-mediated spin-spin entanglement of distant NV 
centers |26| . The ground-state spin triplet can be photoexcited both resonantly (~637nm) and non- 
resonantly to an excited-state spin-triplet orbital doublet. We perform our measurements at cryo- 
genic temperatures (8 K), where these excited-state levels become energetically narrow (27| and their 
fine structureH28l 1291 can be tuned with magnetic, electric and strain fields 11301 I3T1. The NV cen- 
ter's spin is traditionally addressed using its intersystem crossing!; 29j [32| , through which the spin 
is both polarized into the m s = spin sublevel under optical illumination and measured with 
spin-dependent PL intensity emitted in the NV center's phonon sideband (650 — 800 nm). While 
the unique attributes of the NV center's intersystem crossing have made its spin stand out as an 
optically-addressable qubit, the intersystem crossing is not necessary for our optical approaches to 
address the spin. 

We select the subspace spanned by two ground-state spin- triplet sublevels (m s = and +1) as 
our qubit states; the presence of the third sublevel (m s = —1) causes only a small loss in fidelity 
(SI Appendix). We denote these states |0g)and | + l g ). To form the necessary excited state, we 
apply a magnetic field to reach a spin sublevel anticrossing, whose levels are a function of crystal 
strain, spin-spin, spin-orbit, and Zeeman interactions |34j (Methods). The anticrossing we use 
is between the |0 e i) and | + l e \) spin sublevels within the lower-energy excited-state orbital branch 
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Figure 1. A configuration and the NV center A, A configuration within the NV center level 
structure (left), depicting excitation with two optical driving fields from ground states (GS) to excited 
states (ES). At the center of the excited-state anticrossing, the two upper A states \R e i) an d |L e i) 
(bolded green) are the orthogonal, equal superpositions of \0 e i) and | + l e i). An example dark state, 
|D), from the \R e i) A system, is plotted on the rotating-frame Bloch sphere (right), where its polar, 
6, and azimuthal, (p, positions are a function of applied laser power and phase (equation). B, PL from 
resonant excitation as a function of magnetic field and laser frequency illustrating the anticrossing 
between the |0 e i) and | + l e i). C, PL from excitation with two optical fields as a function of the 
detuning of co mw from Scs/h, resonant with either \R e i), \L e \), or centered between both resonances. 

and results in two spin-composite levels (Fig. [Tj3) separated in energy by 8 e \ ~ /z*0.18 GHz. Either 
of these superposed levels, denoted \R e \) and \L e \), can act as the upper state of our A system (Fig. 



In order to address the A-system transitions, we split light from a 637 nm (coi/ (2zr) ~470, 000 GHz) 
laser tunable across the NV center's optical transitions into sidebands (multiples of to mw / (2n) ~4.6 GHz) 
with an electro-optic phase modulator. The relative phase (<p) between the two optical fields that 
are resonant with the A transitions determines the azimuthal position of the dark state on the Bloch 
sphere in the co mw rotating frame. Similarly, the relative amplitude of the two optical fields de- 
termines the dark state's polar angle, 6 (Fig. [l]^). We first observe CPT spectroscopically [1-6] by 
examining the PL under quasi-continuous photoexcitation that optically drives only one of the A 
systems. A sharp dip in PL is observed centered at co mw = Scs/h where Sqs is the mean energy 
splitting between the spin eigenstates (Fig. [lp), indicating that the spin is being coherently trapped 
in the dark state. Because the spin-composite excited states are orthogonal in the |0 e i ) and | + l e \) 
spin subspace, the dark states from each of the separate A systems have opposite azimuthal phases 
but the same polar position on the rotating-frame Bloch sphere for a given optical A-driving con- 
figuration. For this reason, when we tune the laser to equally excite both A systems ("center" 
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curve in Fig. [Tp), their competing dark states quench the PL dip. For subsequent studies, we set 
<*>mw = <5cs/h unless otherwise noted. 

Coherent Population Trapping for Arbitrary-Basis Spin Initialization 
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Figure 2. Time dynamics of coherent population trapping A, Bloch sphere representation of the 
spin state as a function of the CPT interaction time, on resonance with \R e \). Beginning near either 
|0g) (orange) or | + lg) (blue), this process polarizes the spin towards |D) regardless of its initial 
state. Errors are ~3x the point size, and are detailed in the SI. B, Model of the time dynamics using 
a Lindblad master equation approach (description in SI). C, Fidelity of initialized spin state as a 
function of pulse duration. Fidelity is compared to the pure state |D). 



We extend our investigation of the CPT interaction further by probing the time dynamics of the 
resultant spin state. We set the lasers resonant with the \R e i) A system to produce a dark state 
near the Bloch sphere equator. After preparing the initial spin state in either |0g) or | + lg) with 
traditional off-resonant (532 nm laser) optical polarization and microwave electron-spin resonance 
(ESR) techniques [35], we engage the CPT interaction for a variable duration to polarize the spin 
toward the dark state. We then perform quantum state tomography (Methods and SI Appendix) 
on the post-CPT spin state via microwave ESR pulses phase-matched to co mw and subsequent spin 
readout via a second laser resonant with the \0 g ) to (O^) cycling transition! 20] (Fig. [IJ^ and Meth- 
ods). The tomographic reconstructions (Fig. |2]<\) show that the spin state evolves towards the dark 
state regardless of its initial state, and a theoretical model accounting for both A systems (Fig. |2|3) 
is in qualitative agreement with our data (Methods and SI Appendix). As a function of pulse dura- 
tion, the initialization fidelity saturates at about 80% after 100 ns (Fig. |2p)- The fidelity is limited, 
in particular, by decoherence from optical coupling to the second nearby A system, as well as finite 
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T£ spin coherence M36l , spectral diffusion[27j, and some pumping into the third spin state | — l g ) (SI 
Appendix). 




Figure 3. Arbitrary spin-state initialization A, Azimuthal initialization of spins via CPT on reso- 
nance with |R e i). Varying the relative phase between the two optical fields (</>) changes the azimuthal 
location. B, X, Y, and Z projections of azimuthal initialization, on resonance with \L e \j (top) vs. the 
orthogonal state \R e i) (bottom). Error bars are within point size. C, Polar initialization of spins, 
on resonance with \R e i). Varying the relative amplitude between the two optical fields (tan(0/2)) 
changes the polar location. D, Initialization of spins along a great circle canted 7i/ 4 off the polar 
axis, achieved through control of both the relative phase and amplitude of the two optical fields. 
Prior to CPT, the spin was polarized into | + lg) (SI). The CPT pulse duration was 200 ns (A, C, D) 
or 100 ns (B). Errors are ~2x the point size, and are detailed in the SI. 

The allure of this technique is the ability to initialize the spin arbitrarily on the Bloch sphere 
solely by varying the relative phase and amplitude of the two optical fields. In Fig. [3]A we demon- 
strate initialization along different equatorial points of the Bloch sphere by changing the relative 
phase between the two driving optical fields resonant with \R e \). Because \R e i) and \L e i) are or- 
thogonal spin mixtures, tuning the lasers to \L e \) instead is equivalent to shifting cp of the final state 
by tc radians (Fig. [3^3) . Alternatively, by tuning the relative amplitudes of the two optical fields, 
we initialize the spin at various points along a meridian of the Bloch sphere (Fig. |3p)- Finally, we 
combine polar and azimuthal control to demonstrate spin initialization at points along a great circle 
rotated zr/4 radians from the polar axis (Fig. |3p)- 

Arbitrary-Basis Spin Readout via CPT Photoexcitation 

Readout along an arbitrarily-chosen basis ||37| is realized through a complementary process as the 
emitted PL during the CPT interaction with the two optical fields is proportional to the projection 
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Figure 4. Arbitrary-basis spin-state readout A, A configuration recast in terms of ground state 
orthogonal superpositions, the bright \B) and dark |D) states. The driving fields are similarly recast 
as an optical pump on the bright state transition. B, The emitted PL response of the NV center 
spin as it settles into the dark state, starting either near the bright or dark state. This trace is a sum 
of 2.3 x 10 6 iterations with the data binned into 10 ns time intervals. C, Spins initialized at points 
along the equator and read out through DBP. The DBP basis is chosen such that the corresponding 
bright state, indicated in the legend, is at one of four points on the equator (top panel) or one of the 
poles (bottom panel). D, Spins initialized at points along a meridian, mapping out Rabi oscillations, 
and read out via DBP in the same bases as in C. E, (top) All-optical Ramsey experiment, detuned 
such that comw — ^gs/^ = 27r-7.5MHz. The CPT initialization and DBP readout pulses are each 
50 ns in duration, (bottom) Room-temperature Ramsey measurement using ESR pulses with similar 
detuning for comparison. All error bars represent la shot noise. 



of the spin along the bright-state axis. This can be thought of as a recasting of the ground states 
of our A system in terms of the bright and dark states, orthogonal superpositions of the original 
spin eigenstates |0?) and | + lg). The two driving light fields are correspondingly recast as a single 
optical pump acting on the bright state transition (Fig. jS]^), since they do not couple to the dark state 
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from destructive interference of photoexcitation. During the interaction, the spin evolves toward the 
dark state, and the emitted PL provides a measure of the spin state prior to the interaction (Fig. 
|4j3) . This technique, which we refer to as dark/bright-state projection (DBP), bears similarity to 
electromagnetically induced transparency [ 9 1, but we instead measure the transient optical response 
of the NV center rather than the amount of transmitted light. 

To demonstrate arbitrary-basis readout, we prepare the spin state with ESR pulses either along 
various positions on either the equator (Fig. |4p) or a meridian (Fig. [4p) of the Bloch sphere, and 
then use DBP to read out the spin state along six separate bases with bright states corresponding to 
the ±X, ±Y, and ±Z positions of the rotating frame Bloch sphere. The number of photons measured 
is in direct proportion to the projection of the spin state along the chosen axis. The signal-to-noise 
of spin readout using DBP along polar bright states is comparable to traditional spin readout tech- 
niques via the intersystem crossing, while DBP spin readout along equatorial states requires roughly 
3x more averaging (SI Appendix) to achieve a similar signal-to-noise ratio. By combining both CPT 
initialization and DBP readout, we perform an all-optical Ramsey measurement [38 j by varying the 
delay between the CPT and DBP pulses in order to measure the transverse inhomogeneous spin 
coherence time, T| (Fig- |4p!, top). Collapses and revivals in the signal are indicative of hyperfine 
coupling to the 14 N spin. The all-optical response is similar to Ramsey measurements taken at room 
temperature using ESR pulses and traditional intersystem crossing-based initialization and readout 
(Fig. |4p;, bottom). 

Arbitrary- Axis Spin Rotations via Stimulated Raman Transitions 

Within this same optical coupling framework, we also demonstrate unitary spin rotations about 
any qubit axis via SRT By detuning coi from resonance while keeping co mw = 5Qs/h, driving 
the A system produces adiabatic energy shifts of the bright state during the laser pulse without 
modifying the dark state energy, generating unitary spin rotations 171 l8l IT5l [191 along the dark/bright 
state Bloch sphere axis. In order to drive rotations about an equatorial axis, we tune the two equal- 
intensity (tan(0/2) = 1) driving fields to be centered between the \R e i) and \L e \) resonances, such 
that SRT generated from both A systems add constructively while CPT effects from both reduce 
coherence but produce no net spin polarization due to competing dark states. In Fig. [5^, we 
present the dynamics of SRT spin rotations along two different equatorial rotation axes of the Bloch 
sphere {"ax" or "cry"), corresponding to different relative phases (</>) of the two optical fields. We 
measure a 69% process fidelity for a "ax" or "ay" 7r-rotation, limited largely by spontaneous decay. 
Rotations about non-equatorial axes, such as the polar axis l|T9l ("az"), are also achievable in this 
system (Fig. |5|3) but require different configurations of the light fields (SI Appendix). 

Finally, to illustrate the full suite of these optical control protocols, we present an all-optical Hahn 
echo measurement of an NV center spin's homogeneous spin coherence time, T2. This measurement 
consists of a CPT laser pulse for spin initialization along the Bloch equator, followed by a SRT laser 
pulse to flip the spin to produce an echo, and finally a DBP readout pulse to measure the final 
spin state along an equatorial basis (SI Appendix). We determine T2 ~900 us, corroborated by an 
ESR-based Hahn echo measurement at room temperature (Fig. |5p). 
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Figure 5. All-optical control of the NV center spin A, Bloch sphere representation of "ax" and "cry" 

coherent rotations at ~10MHz due to SRT. The two measurements correspond to different relative 
EOM driving phases (</>), separated by Tt/2 radians, and show the trajectory of a spin originating 
near |0g). (orange) and | + l g ). (blue). The axes of rotation are added as guides to the eye. B, Bloch 
sphere representation of "cr z " coherent rotations due to SRT, showing spin trajectories originating 
near orthogonal points on the equator (maroon and grey). For A. & B., errors are ~2x the point size, 
and are detailed in the SI. C, All-optical Hahn echo measurement (green points) consisting of CPT 
spin initialization, a re spin rotation via SRT, and DBP spin readout fit according to the equation on 
the graph. Room temperature Hahn measurement via ESR pulses is also shown (grey points). Error 
bars are lcr shot noise. 



Conclusions 

We demonstrate all-optical initialization, readout, and coherent unitary rotations of an individual 
NV-center spin, forming a triumvirate of protocols for single-spin control that can be performed 
along any arbitrarily-chosen basis. Using these protocols, we demonstrate two measurements of 
transverse spin coherence solely with optical pulses. The ability to select any basis allows for quan- 
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turn operations to be implemented directly without the need for extra control steps to project onto or 
from the preferred energy eigenstate basis. This eliminates the need for ESR operations ||35| , enabling 
control of individual spins within a much smaller device footprint, with promise for large-scale im- 
plementations of spin arrays|39| or photonic networks[21, 22]. Perhaps most importantly these 
methodologies mitigate the need for the NV center's intersystem crossing spin-selectivity and thus 
can be used to investigate and control a wide array of defects and other localized quantum states 
in solid-state materials, not just those with NV-like structures BUI |41). As such, these techniques 
open the door to exploring quantum coherence and developing quantum information platforms in 
a broad range of semiconductors and nanostructures. 

Methods 

Sample - The sample was a 2 x 2 x 0.5 mm electronic grade diamond purchased from Element 
Six, consisting of < 5 ppb nitrogen that was irradiated with a lel4 electrons /cm 2 , 2 MeV dose 
and subsequently annealed at 850 °C for two hours. Ti/Pt/Au devices, consisting of DC pads 
and a short-terminated waveguide, were deposited on the sample using standard photolithographic 
techniques. All experiments were performed in a confocal microscopy setup (SI Appendix) with 
a liquid helium flow cryostat held at 8 K. The sample was thermally sunk to the cryostat and the 
waveguide was wirebonded to a microwave line in the cryostat for on-chip ESR. The studied NV 
center excited state orbital strain splitting between m s = spin sublevels varied from 4.6 GHz to 
5.8 GHz between cryostat cooldowns as thermal cycling modified the crystal strain. As a result, 
the DC-applied magnetic field at which the lower-branch excited state spin anticrossing occurred 
varied (550 — 750 G) which led to variations in the ground state spin splitting between \0 g ) and 
| + l g ) (S GS /h ~4.3 to 5.0 GHz). 

Experimental Techniques - Our confocal microscopy setup consists of a 100 mW 532 nm non- 
resonant excitation laser used for the standard NV center intersystem crossing spin initialization 
and readout protocol lfT7) , and two tunable 637 nm lasers resonant with various NV center optical 
transitions. The light field from one of the resonant lasers was fiber-coupled to an EOM in order to 
split the optical field, at coi, into different frequency sidebands, separated by w mw , to optically drive 
the A system. For the CPT and SRT measurements, a second resonant laser functioned as a one-color 
spin-state readout laser along the Sz basis by being resonant with the \0 g ) to \0 e 2)- transition! 20 1 
(Fig. [IJ^) resulting in higher collected PL when the spin was in |0g). In Fig. |5j: only, the light field 
from this second resonant laser was instead fiber-coupled to a second EOM, where the first laser 
was used to perform CPT and DBP and the second laser was used to perform SRT for this Hahn 
echo pulse sequence. All three lasers were gated using separate acousto-optic modulators (AOMs) 
for pulse timing control. They were subsequently passed through a variety of polarization optics, 
combined with beamsplitters, and focused onto the sample with a 0.85 numerical aperture 100 x 
microscope objective that is aberration-corrected for the cryostat window. PL from the NV's red- 
shifted phonon sideband was collected back through the objective, filtered by dichroic beamsplitters, 
and focused onto a silicon avalanche photodiode (APD). 

Microwaves to drive the EOM(s) and for on-chip microwave ESR driving [ 35 1 originated from the 
same signal generator at frequency co mw / (27r), which varied from 4.3 — 5.0 GHz due to variations in 
Sqs I h from changes in NV center strain. The microwaves going to the EOMs and to the sample for 
ESR passed through IQ modulators for phase control between the various CPT, SRT, DBP, and ESR 
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pulses. These microwave signals were also gated in time and amplitude-controlled using microwave 
modulators and switches. Timing for the microwave switches, AOMs, and IQ modulators were 
controlled by an arbitrary waveform generator, a PulseBlaster card, and a pulse-pattern generator. 
Pulse sequences used for these experiments consisted of a traditional initialization pulse at 532 nm 
as a spin reset, followed by a sequence consisting of a number of the following techniques: on- 
chip microwave ESR pulses as well as techniques utilizing the 637 nm tunable lasers, including CPT 
spin-state initialization, DBP spin-state readout, \0&) spin-state readout, and/or SRT coherent spin 
rotation. Details for each pulse sequence are discussed in the SI Appendix. A magnetic field was 
applied along the NV center axis with a permanent magnet on a motorized stage and was adjusted 
to tune to the anticrossing used. 

Quantum State Tomography - To perform quantum state tomography on our CPT spin-state ini- 
tialization and SRT coherent rotation, we read out the X, Y, and Z projections of the post-interaction 
state. All projections were mapped onto the Sz basis using ESR pulses and then read out with the 
laser resonant with \0 e 2). We applied a Bayesian approach to the tomographical reconstruction of 
the spin state [42 j, detailed in the SI Appendix, that takes into account finite readout contrast, laser 
drift, and axial/length imperfections in the microwave rotations used to project the different spin 
components. 

Theoretical Modeling - To describe the dynamics of the NV center spin under optical excitation in 
the A level configuration [ 43 1, we include five energy levels: two out of the three ground-state levels 
|0g), | + l g ), the two mixed excited states \L e i) and \R e i) as the upper state of each A system, as 
well as the intermediate singlet |S) which here plays a role for unintentional intersystem crossings. 
The Hamiltonian, in the rotating frame, for the subspace spanned by these five basis states can be 
expressed as 

H = X> a |«)<«|+ £ E (pGE\E A )(G g \+h.c.) t (1) 

a G=0,1E=R,L 

where the first sum runs over all states a = g , +l g , L e \, R e \, S with corresponding energies eo g = e+i ? 
(where is detuning of coi from resonance to a A system), €R el = 0, ei el = —5 e \, and e$- The laser 
excitation from one of the lower states G = 0, 1 to one of the upper states E = L, R is described by 
the Rabi frequencies in the rotating frame, 

Qie = Qcos(0/2) (2) 
n 0E = ±Qsin(0/2)e^ (3) 

where the upper (lower) sign holds for E = R (E = L). 

We studied the time evolution of the system by numerically solving the Lindblad master equation 
45 1 for the density matrix of the NV center in the rotating frame. In addition to coherent processes 
such as excitation from the two driving fields, the master equation also accounts for spontaneous 
decays of charge and spin with some rates known from independent experiments. In the idealized, 
long-time limit case, with only one excited level included, the resulting eigenvector with eigenvalue 
corresponds to the dark state: 

\D) =cos(0/2)|Os)-e T, '*sin(0/2)| + l g ) (4) 
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where the upper (lower) sign holds for the single excited state level being E = R(E = L). In 
actuality, the steady state is described by a mixed state which can deviate slightly from \D){D\. The 
simulated behavior of the NV spin during CPT and SRT is in good qualitative agreement with the 
experimental data (Fig. |2]and Fig. S6 S9b. Further details can be found in the SI Appendix. 
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Experimental details 

A schematic of our confocal setup is provided in Fig. [Sllthat incorporates a continuous wave (CW) 
green diode laser (532 nm), two CW tuneable red (637 ran) diode lasers with optional sideband 
wavelengths generated via electro-optic phase modulators (EOMs), NV center photoluminescence 
(PL) collection (650 — 800 nm) via an avalanche photodiode, and on-chip microwave electron spin 
resonance techniques (ESR) all gatable in time. All timing sequences used were programmed into 
an arbitrary waveform generator and are described later in this Supplementary Information. 

While the effects we presented in the main text were observed in multiple NV centers, care was 
taken to select an NV center that had several desirable properties. We chose an NV center with 
a reasonable T| time (~1 |a.s) due to the finite duration of the CPT/SRT/DBP pulses. Because the 
energy separation between \R e i) and \L e \) A systems was quite small {5 e \/h — 0.18 GHz), it was 
desirable to choose an NV center with a narrow inhomogeneously-broadened optical linewidth 
(0.05 GHz FWHM for the NV center presented). We also chose an NV center that had a transverse 
strain splitting which varied between 4.6 GHz and 5.8 GHz as measured between the two orbital 
m s — spin sub levels. The variations in this NV center's strain occurred between cooldowns due 
to thermal cycling of the cryostat. The strain was in a range which allowed us to produce an 
avoided level crossing (anticrossing) in the lower excited state orbital branch between the m s — 
and m s — +1 spin sublevels with the application of an external magnetic field along the N-V axis, 
(550 — 750 Gauss). The resulting ground state spin splitting at these magnetic fields was relatively 
high (Sqs /h =4.3 — 5 GHz) allowing for the EOM sidebands to be widely spaced, meaning the 
unused laser harmonics generated by the EOM were far from any NV center resonances. 
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Figure SI. Schematic of experimental setup. Diagram detailing the optical excitation paths and 
photon collection path along with the microwave electronics, timing electronics, cryostat, and mag- 
netic field as described in Methods. A PL scan of the region 6 |a.m below diamond surface, where the 
NV centre investigated is located (within the smoke-colored triangle). A short-terminated on-chip 
waveguide wire used to apply microwave ESR pulses for ground state spin manipulation is visible 
in the lower left. Deposited metallic pads on the right and top of the image are for applying dc 
voltages to the sample to affect the orbital splitting of NV center if necessary, but were not used in 
the present experiment. 



It is also possible to generate an anticrossing in the upper excited state orbital branch but the re- 
sulting anticrossed eigenstates are much closer together in energy, making it prohibitively difficult 
to couple to an individual resonance. Off-axis magnetic fields could be used to increase the split- 
ting between anticrossed levels, but they would also reduce the spin-selectivity of the intersystem 
crossing (ISC) used for comparison to the CPT-based readout. 

In order to fully control the phase and amplitude of the CPT laser system, we adjust properties 
of the microwaves driving the EOM. To adjust the dynamic phase between the two used light fields, 
we control the phase of the co mw microwaves driving the EOM with an IQ modulator, which moves 
the dark state azimuthally about the rotating-frame Bloch sphere. We adjust the relative amplitude 
of the two used laser driving fields by varying the microwave power driving the EOM using a 
mixer, allowing us to adjust the polar position of the dark state on the Bloch sphere. While the 
relative amplitude of the two resonant optical fields adjusts as expected, the summed amplitude of 
the two A-resonant sidebands also changes. This is because the EOM splits light at a>i into several 
harmonics (of which we only select two) separated by co mw , with amplitudes of the harmonics being 
a function of the amplitude of EOM microwave driving power in the form of a Bessel function. To 
correct for this, we use an optical variable attenuator to compensate for this overall amplitude 
variance to within an order of magnitude in an attempt to keep Q (Fig. 1A) fixed. In addition, 
the overall frequency of the tunable laser, a>i, would drift, and due to the small separation in the 
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anticrossing between \R e i) and \ L e i), S e i/h — 0.18 GHz, we needed to occasionally recenter the laser 
frequency to the appropriate tuning on the order of every 10 minutes with 0.02 GHz laser frequency 
resolution. 



Quantum state tomography of arbitrary initialization and rotation 



In order to analyze our initialization and control protocols, we performed Bayesian quantum state 
tomography to characterize the various process output states and compute the corresponding fi- 
delities found in the main text. This approach allows for an accurate statistical inversion of repeated 
projective measurements that are subject to both stochastic and systematic error. In contrast to max- 
imum likelihood estimation, this approach always yields both point estimates and corresponding 
error bars that are physical for states near the boundaries of the allowed state space; moreover, it 
relaxes the assumption of asymptotic normality, achieving consistent estimates in the face of a finite 
number of measurements [42|. 
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Figure S2. Pulse sequence for arbitrary initialization and rotation. The above pulse sequence was 
used for the data presented in Main Text Fig. 2, Fig. 3, and Fig. 5A,B. All experiments investigating 
these protocols consisted of ~ 10 6 iterations of this pulse sequence to achieve a sufficient signal-to- 
noise ratio. 



To individually study CPT or SRT, the NV center spin state is first prepared using non-resonant 
532 nm laser light to both mitigate photoionization and polarize the NV center ground state spin 
into m s = via ISC decay with roughly 75-80% fidelity [32]. |46|. The state is then prepared on 
various places of the Bloch sphere with an ESR "preparation pulse" before either a CPT or SRT red 
laser pulse polarizes or rotates the spin, respectively. After the red laser, an additional ESR "QST" 
pulse is applied to rotate the X, Y, or Z spin projection onto the |0g)/| + lg) readout measurement 
basis, phase-synced with the microwaves driving the EOM. The timing of this QST pulse is chosen 
to coincide with the constructive rephasing of the three 14 N hyperfine Larmor frequencies, cor- 
responding to a delay of 450 ns. Alternatively, a 7r-pulse can be added before the QST pulse to 
induce an echo of the spin coherence rather than use this rephasing. Finally, the spin state is read 
out along this basis by measuring PL intensity during a single-color red laser resonant with the 
|0g)/|0 e 2) cycling optical transition 1521 . 

As a normalization, the state is more fully initialized into |0g) (or after 532 nm excitation 

by subjecting the spin to 637 nm laser light for 1 us resonant with the opposite | + (or \0g)) spin 
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state and one of the A systems. This scheme depletes the aforementioned optically-driven sublevel 
in the ground state and populates its counterpart \0g) (or The readout contrast, C, between 

the |0g) and \+lg) sublevels increases from about 62% from ISC spin polarization alone to about 
84% after this purification. Because none of the CPT (or SRT) interactions appeared in practice to 
polarize the state beyond this level, this purified contrast serves as a consistent normalization for 
each state reconstruction and fidelity. 

To begin our analysis, in terms of the density matrix, p, the expectations are 

(X) = Tr ( a z U Y n pU f Y n 

<Y>=Tr(cr z L7x 5 pU x? ) 
(Z)=Tr(cr 2 p), 

while the expected fluorescence levels (in photon counts) are defined by the resonant laser normal- 
izations, 
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(^{Y}> 


= ^i) 


^<z>> 





We treat each of the data, D k , as subject to normal error a k from the model prediction (F) k 
whose other parameters (J 7 ^), C, and those described below) are contained in a vector X such that 

prob(D|^,X) =n-^exp f~ (Pfc ijf^ l ■ (SI) 
k V27tcr k \ 2a{ J 

Because of shot noise, the lower limit to a k is about 1/ (F) k but due to both the drift of the laser and 
stage mechanics, cr k was about a factor of two to five larger in practice. To capture this uncertainty 
in the expected mismatch, we set 



prob (a k \a k ) = -^exp [ --\ ) , (S2) 
V not 




so that we marginalize each cr k around a region of order W k , expressing that a k should be on the 
order of, but not necessarily equal to, a k ; this sort of construction makes our estimation statistically 
robust against data of unusually large drift (outliers). After performing the integration over all 
positive a k , we have that 



prob (DI^X) = ( V2n<r k ( 1 + Jf )fc 



24 




(S3) 



In our analysis, we set c k = 1\J (F) k . This likelihood is also used for normalization parameters ^\Q g ) 
and C in the posterior probability density to infer them simultaneously with the axial projections 
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and systematic errors. 

In addition to the random noise of the experiment, the microwave pulses used to rotate the X 
and Y components to the Z axis for readout suffer from small systematic errors in their relative 
phase, which creates an offset from the proper rotation axis, and also in their duration, which 
creates an offset from the proper rotation length. The unitary operators above thus deviate from 
the ideal case and can be redefined, following Dobrovitski et al. [47J, by 
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and a is a vector com- 



where n x = (l,e y ,e z ) / (l + e 2 . + e 

posed of the Pauli matrices. Here, the e and v terms are the axial offsets and cp and 9 are the length 
offsets, all of which are measured in units of angle. From calibration of the IQ modulator and the 
discrete nature of the delay generator that governs the length of the pulses, we estimate that these 
factors are at most 5° in angular error and somewhat conservatively set the prior densities for each 
of these terms as a normal density of mean zero and standard deviation of 5°. The effect of this 
correction is most noticeable from the fact that the X and Y axis projection estimates almost always 
have larger uncertainties than the Z axis estimate. 

Lastly, we use the non-informative reference prior[48J for the density matrix, 



prob (p) = 0.00513299 (1 - r 



-1/2 



log 



L(l + r) 



sin 9, 



(S4) 



where r and 9 are the standard spherical coordinates for the Bloch vector used to parameterize p. 
As expected, given the large quantity of data collected, the choice of prior had no discernable effect 
on our inferences. 

To obtain the marginal densities used for the point estimates and error bars of the projections 
and fidelities found in the main text, we use MT— DREAMzg/ a Markov Chain Monte Carlo tech- 
nique, to sample from the 11-dimensional posterior probability density [49|. This technique uses 
multiple random walk chains of a multiple-try Metropolis-Hastings rule |50l applied to an adaptive 
proposal distribution generated from past samples to generate new samples in accordance with the 
posterior probability density. We found the chains to have fast convergence and good mixing prop- 
erties, and the corresponding sampler output to have low autocorrelation using four independent 
chains and a multiple-try parameter of seven. The fidelities for experimental data are obtained by 
calculating the fidelity between the Bloch vectors from the random walk sampler and an ideal vector 
of unit length pointing along the same axis as the corresponding random walk sample. In this way, 
the fidelity is computed with respect to a perfect rotation of the state generated from the resonant 
laser pumping scheme described above. The point estimates of the projections and fidelities are 
the mean of the respective marginal densities while the error bars are the highest posterior density 
68.2% credible intervals ITST1. 



In Figs. S3 S4 and |S5j we present the corresponding X, Y, and Z projections for the rotating- 
frame Bloch spheres describing arbitrary initialization via CPT presented in the main text, as well as 
some additional sets of projections and Bloch spheres not presented in the main text. As described 
in the main text, all of these qubit Bloch spheres are within the co mw = 5q$ I h rotating frame. This 
is because our standard qubit states |0g) and are separated in energy by 5qs, and so each 

spin state precesses at the Larmor precession frequency of Sqs/Ti. By reading out the X, Y, and Z 
projections with ESR pulses phase-matched relative to the same co mw — 5q$ / fi angular frequency, 
we capture a snapshot of a fixed point on this rotating-frame sphere. Conversely, in a system where 
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the qubit states are not energy split, the Bloch sphere would not be precessing. A CPT initialization 
laser pulse length of 200 ns is used for Figs. S3j S4j and S5 In some instances, noted below, this is 
not long enough to fully polarize the spin. 
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Figure S3. Projections for azimuthal initialization of spins. The azimuthal location of the polar- 
ized spin state is rotated along the equatorial plane by varying the relative phase between the two 
colors. X, Y, and Z projections are plotted on the left, and reconstructed Bloch spheres are plotted 
on the right. Top: Prior to the CPT interaction, the state was |0 g ). Bottom: Prior to the CPT inter- 
action, the state was in | + l g ); this data is presented in Fig. 3A in Bloch sphere form. Error bars on 
projections are the 68.2% highest posterior density credible intervals from the Bayesian analysis. 



Note that in Fig. S4 there is an asymmetry between the Z projection of the dark state near the 



poles depending on whether the state was prepared in |0g) or This is partly due to the fact 

that a longer initialization pulse would be required to fully move the state to the opposite pole. In 
addition, some of this imbalance could also be due to differential spin coupling of the ISC and its 
resultant decay before readout occurs. This effect can also be seen in the Fig. S5 data, where we 
combine phase and amplitude control of the light fields in order to initialize points along a great 
circle rotated n/ 4 off the equator. 

Finally, we should also note that the collection of Bloch sphere representations in all Main Text 
and Supplementary figures are not all viewed from the same vantage point, some were rotated to 
better show the results. However, within a single figure the angle of view is fixed. One angle of 
view is used for the Bloch spheres in Fig. 2 and Fig. |S6| (Time dynamics of CPT initialization). A 
second angle of view is used for the Bloch spheres in Fig. 3, Fig. S3 S4 and |S5] (Arbitrary spin-state 
initialization). A third angle of view is used for the Bloch spheres in Fig. 5, Fig. |S7j |S8j and S9 
(All-optical control of the NV center spin). 
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Figure S4. Projections for polar initialization of spins. We vary the polar location of the polarized 
spin state by varying the relative amplitude between the two colors. X, Y, and Z projections are 
plotted on the left, while the same points in a Bloch sphere representation are plotted on the right. 
Top: Prior to the CPT interaction, the state was in |0 g ). Bottom: Prior to the CPT interaction, the 
state was in | + l g ); these data are also presented in Fig. 2B in Bloch sphere form. Error bars on 
projections are the 68.2% highest posterior density credible intervals from the Bayesian analysis. 
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Figure S5. Projections for initialization of spins along an off-axis great circle. Here we vary both 
the relative amplitude and phase between the two colors to place the spins at points along a great 
circle, tilted jt/4 off of the equator. Top: Prior to the CPT interaction, the state was in |0 g ). Bottom: 
Prior to the CPT interaction, the state was in | + l g ), this data is presented in Fig. 3D. Error bars on 
projections are the 68.2% highest posterior density credible intervals from the Bayesian analysis. 
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Model of arbitrary initialization and rotation 



The Hamiltonian describing our system (Equation [1] from Methods) is presented below in matrix 
form, 



H = 
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where the ordering of the states in the matrix is: {| + lg), \0g), \R e \), \L e i), \S)}, Ai is the detuning 
of the laser frequency (col) from resonance to the \R e \) A system, 5 e \ is the separation of the excited 
state levels, Q is the optical Rabi frequency, (p is the relative phase between the two coherent light 
fields, and tan(0/2) is the relative amplitude between the driving fields. As such, <p and 8 will 
describe the azimuthal and polar angle, respectively, of the resultant dark state. 

The time evolution of the system is described by the Lindblad master equation |44l |45| , 



i[p,H] 



0L,tx! 
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Wp, 



(S5) 



with the Lindblad operators a aa — |a)(«| = <v a <V a and cr a / a = cr , = |a')(a|. For n — 5 levels, the 
density matrix p is a Hermitian 5x5 matrix and can thus be described by n 2 — 25 real parameters 
(n 2 — 1 = 24 including the normalization condition Tr (p) — 1). The superoperator W can thus be 
viewed as a 25x25 matrix with rank 24. We denote the decay rate from the excited states (E = L,R) 
to the ground states (G = 0, 1) with T — T£ eli Q g , the rate for inter-system crossing from the excited 
states to the singlet Tj = Te Ai s, and the inverse intersystem crossing rate from \S) to one of the 
ground state levels as T- = Ts,g*- The spin relaxation rate in the ground state is T\ = X/T\ = r+i„,o x 
and at sufficiently low temperature Io +i ~ 0. The pure dephasing between the two ground state 
levels is denoted 7 = I/T2 = Tq q . All other rates are set to zero. 

The state of the system after optical excitation during time t is obtained as 

Wf, 



p(t)=e™p(Q), 



(S6) 



where we choose one of the ground states as the initial state, p(0) — [Og) (0^| or p(0) = | + lg) (+1 



We typically determine p(t) by performing the exponentiation Eq. (S6| numerically. The resulting 
Bloch vector in the ground state subspace can be obtained from 



b(0=Tr(<rp(0), 

where the components of tr are the Pauli matrices in the ground-state subspace, 



a x = 
a y = 
a z = 



\+lg)(0g\ + \0g)(+lg\, 

i(\ + h)(0g\-\0g)(+h\)' 
\0g)(0 g \-\ + lg)(+lg\. 



(S7) 



(S8) 
(S9) 
(S10) 



In the idealized case Ti = 7 = Tj = 0, and with only one of the excited levels included, the 
stationary state p in the long-time limit t > 1/F obtained from p — as the eigenvector of W with 
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eigenvalue is the dark state: 

\D) =cos(0/2)|O g ) -exp(Ti0)sin(0/2)|+l g ) (Sll) 

where the upper (lower) sign holds for the single excited state level being E — R (E — L). The 
fidelity of the state after a finite pumping time t with realistic parameters is then 

F(t) = (D\p(t)\D). (S12) 

The experimentally obtained fidelity is shown in Fig. 2C in the main text. 

We use this model to simulate the time evolution of the Bloch vector b(f) during 500 ns of CPT 
initialization plotted in Fig. 2B in the main text (Fig. [S6] l and the 200 ns of SRT rotation. We fix the 
excited state splitting b e \ — 180 MHz, and T' using the~known lower singlet lifetime of 371ns|32|. 
We also fix the rate for both the \R e i) and \L e \) states into the singlet, T,, as one half of the known 
ISC rate of m s — ±1 spin states to the singlet of ~ 74 MHz, since \R e i) and \L e \) are composed 
of equal mixtures of m s = and m s = +1 states. The detuning of the optical fields relative to 
\R c i), Ai, is fixed to the experimentally measured detuning for each operation (CPT initialization, 
SRT a x rotation, SRT cr y rotation, SRT <7 Z rotation). Using a weighted least squares approach, we fit 
the overall driving amplitude, Q, relative amplitude, tan(0/2), relative phase of the driving fields, 
(p, and the decay rates T, Y\, and 7, described above. All simulations show qualitative agreement 
with the experiment, however certain traces appear to have out-of-phase behavior of individual 
projections, marginal agreement with the initial state, or other disagreements. It appears that the 
model captures the essential physics but cannot fully account for certain ill-defined nuances such 
as the transients during the turn-on/ off of the optical fields and effects related to the hyperfine 
spectrum and T£. Therefore, the fitted values for the decay parameters are skewed by effects not 
considered in the model. A full set of fixed and fit parameters is found in Table [Slj 

As mentioned in the main text, it may be possible to rotate about any arbitrary axis, but there 
are a few considerations to be made. SRT are a dispersive interaction whose strength is propor- 
tional to 1 / A^, whereas CPT is an absorptive process whose interaction strength is proportional 
to 1/A?. To take advantage of SRT, sufficient detuning is necessary to diminish absorptive effects 
that are non-unitary and cause the spin to polarize along the dark state rather than rotate. In the 
case of rotations about an equatorial axis, we take advantage of the two competing A systems by 
tuning oj 1 exactly between \R e i) and \L e \) resonances. These A systems have opposite equatorial 
bright states and we detune from both in opposite directions, essentially causing SRT effects to add 
constructively and the CPT effects to add destructively. In order to rotate about an axis off of the 
equator, care must be taken to couple more strongly to one of the A systems than the other because 
their corresponding bright states are no longer orthogonal, as the bright states for each A system 
will have orthogonal azimuthal phases but the same polar component. As discussed above, a more 
widely spaced anticrossing to decrease the competition between the two upper states would aid in 
ensuring unitary non-equatorial axial rotations. The special case of rotation about the polar axis, 
which we demonstrate in Fig. 5B, was also shown in Buckley et fl/. |T9| but not in the context of a 
A system. In our case, we sufficiently detune ~450 MHz from the \R e i) such that any effects from 
CPT are greatly diminished. 
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CPT initialization: Tomographic Projections 
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Figure S6. Time dynamics of arbitrary initialization: theory vs. experiment. As a function of CPT 
initialization pulse duration, tomographic reconstructions of the spin state are plotted alongside a 
simulation of the resultant state using the model. Top: Prior to the CPT interaction, the state was 
|0 g ). Bottom: Prior to the CPT interaction, the state was in | + l g ). Error bars on projections are 
the 68.2% highest posterior density credible intervals from the Bayesian analysis. Parameters to 



simulate spin-state initialization are found in Table SI 
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Figure S7. Time dynamics of arbitrary coherent cr x rotation: theory vs. experiment. As a function 
of the duration of a SRT ax rotation pulse duration, projections of the resultant spin state are plotted 
alongside a simulation of the resultant state using the model. Top: Prior to the SRT interaction, the 
state was |0 g ). Bottom: Prior to the SRT interaction, the state was in | + l g ). Error bars on projections 
are the 68.2% highest posterior density credible intervals from the Bayesian analysis. Parameters to 
simulate the ox rotation are found in Table |S1 
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Figure S8. Time dynamics of arbitrary coherent cry rotation: theory vs. experiment. As a function 
of the duration of a SRT oy rotation pulse duration, projections of the resultant spin state are plotted 
alongside a simulation of the resultant state using the model. Top: Prior to the SRT interaction, the 
state was |0 g ). Bottom: Prior to the SRT interaction, the state was in | + l g ). Error bars on projections 
are the 68.2% highest posterior density credible intervals from the Bayesian analysis. Parameters to 
simulate the cry rotation are found in Table [SI 
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"cr " rotation: Tomographic Projections 
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Figure S9. Time dynamics of arbitrary coherent t7 Z rotation: theory vs. experiment. As a function 
of the duration of a SRT cz rotation pulse duration, projections of the resultant spin state are plotted 
alongside a simulation of the resultant state using the model. Top: Prior to the SRT interaction, 
the state was |X g ). Bottom: Prior to the SRT interaction, the state was in | — X g ). Error bars on 
projections are the 68.2% highest posterior density credible intervals from the Bayesian analysis. 



Parameters to simulate the c z rotation are found in Table SI 
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Parameter 


CPT initialization 






fz 


Fixed Level 


&e\ 


180 MHz 


180 MHz 


180 MHz 


180 MHz 


Parameters 


A L 


-0.684 MHz 


-90 MHz 


-90 MHz 


-450 MHz 


Optical 


Q 


46.507 MHz 


62.021 MHz 


62.756 MHz 


84.104 MHz 


Driving Field 


e 


1.708 rad 


4.774 rad 


1.763 rad 


n rad 


Parameters 


<P 


0.395 rad 


4.152 rad 


2.683 rad 


0.424 rad 


Decay 


r 


35.114 MHz 


17.115 MHz 


19.719 MHz 





Parameters 




37 MHz 


37 MHz 


37 MHz 


37 MHz 




r; 


2.701 MHz 


2.701 MHz 


2.701 MHz 


2.701 MHz 






0.373 MHz 













7 











28.459 MHz 


Initial State 


ta 


0.640 


0.839 


0.852 


0.602 


Parameters 


e A 


0.164 rad 


0.327 rad 


0.347 rad 


1.844 rad 






2.526 rad 


4.705 rad 


2.850 rad 


1.471 rad 




?b 


0.649 


0.977 


0.752 


0.621 




e B 


2.788 rad 


2.325 rad 


2.774 rad 


1.870 rad 




<pB 


3.122 rad 


3.450 rad 


3.472 rad 


4.425 rad 



Table SI. Simulation Parameters. Simulation parameters for CPT initialization, C7 X/ cy, and <7 Z . In 
the case of CPT initialization, C7 X/ and cy, the initial state parameters state A refers to |0 g ) and state 
B refers to | + l g ). For a z , state A corresponds to |X g ) and state B corresponds to |— X g ). 



S15 



Sources of decoherence 



The closeness in energy of the two anticrossed eigenstates, |jR b i) an d \L e l), causes a loss in fidelity. 
This is because even when resonantly tuned to a single eigenstate, there is still off-resonant coupling 
to the other A system. As the phases of the resultant dark states between the two are orthogonal, the 
overall length of the vector pointing to the final state within the Bloch sphere is reduced. The model 
shows good qualitative agreement to our presented data, revealing a similarly mixed state as a result 
of the competing A systems. Therefore, this competition is one of the most significant source of 
decoherence in our measurement. Finding a more widely spaced anticrossing would help alleviate 
this issue by decreasing the coupling to the other state; within the model, a purer final state (higher 
fidelity) results with increasing separation of the two excited states. Experimentally, one method 
to achieving more widely spaced anticrossing would be to slightly misalign the field because off- 
axis fields will increase the separation in any anticrossing. However, it should be noted that a 
misalignment of the field will also change the eigenstates, possibly producing more transitions 
from spin-mixing into the |— l g ) spin sublevel, which also causes a loss in fidelity. 

Further sources of decoherence include the transverse inhomogeneous spin coherence time 
T| ~1 [is as well as the fact that our spin sublevels are further split into three nuclear hyperfine 
states due to the 14 N in our NV center, with each transition split by ~2 MHz. Since our measure- 
ments are averages over ~ 10 6 replications of the same experiment, any individual replication has 
an equal probability of being in any of the three hyperfine states, causing our selected co mw to be 
~ ±2 MHz detuned from the actual ground state splitting two thirds of the time. We note that no 
nuclear polarization was observed at this excited state anticrossing at cryogenic temperatures. This 
hyperfine spectrum and finite T| effectively set limits on how long our A interaction remains phase 
coherent with the spin, reducing transverse coherence of the interaction. The resultant steady state 
is a balance between the strength of the interaction and these decoherence mechanisms (along with 
the other mechanisms mentioned in this section). As such, an echo sequence would not eliminate 
these effects while the interaction is taking place, but only during the rest of the measurement se- 
quence. In addition, as the NV center is a solid-state defect, spectral diffusion causes a broadening 
of the natural linewidth of the resonances |27], which will also contribute to a lower fidelity. 

Finally, since we are examining a A system contained within a more complex level structure, 
the spin will end up outside of our qubit subspace a fraction of the time. First of all, the spin can 
end up passing through the long-lived ISC, which is accounted for in the model as a singlet level, 
and does manifest as a decoherence mechanism. Secondly, since our qubit states, |0g) and | + lg), 
are a subspace of a spin-triplet ground state, the spin can transition to | — lg) some of the time. This 
is because, experimentally, the excited state levels \R e i) and \L e \) each contain anywhere from 1% 
to 3% of the excited state spin sublevel | — \ e \). A spin in the|— lg) sublevel is not resonant with 
our red lasers and will therefore appear dark. This suggests that there might be higher fidelities 
resulting from implementation of this type of control within a system that lacks decay mechanisms 
such an intersystem crossing. 

Through our tomographic reconstructions, we are able to determine that the fidelity of our spin- 
state initialization saturates at roughly 80% after 100 ns (Fig. S10A ). The fidelities of a a\ or n 
rotation are as high as 69% if we assume a pure initial state as the initial state instead the mixed 
state used in the experiment (Fig. SlOp , C). A c z n rotation could have process fidelities as high as 
90% (Fig. Slop ), but as the initial states used for that experiment were of lower fidelity than any 
other experiment, this estimation may be rather optimistic. 
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Figure S10. Fidelities of initialization and rotation. A, Fidelity of initialization as a function of 
pulse duration for an initial state |0 g ) and | + l g ). B, Process fidelity of o~x rotation compared to a 
perfect n rotation, as a function of pulse duration for an initial state |0 g ) and | + l g ). C, Process 
fidelity of cry rotation as compared to a perfect n rotation, as a function of pulse duration for an 
initial state |0 g ) and | + l g ). D, Process fidelity of cr z rotation as compared to a perfect n rotation, 
as a function of pulse duration for an initial state |X g ) and |— X g ). For all process fidelities, we 
compare the resultant state, as a function of pulse duration, to a state exactly n out-of-phase with 
the initial state about the rotation axis. We renormalize the initial mixed state to a pure state in 
order to compute the fidelity loss from the CPT or SRT process alone. 
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Arbitrary spin-state readout 



The following pulse sequence was used to perform the DBP readout protocol. 



Traditional 
Initialization 



532 nm Laser 



DBP Readout 



DBP 



Tunable 637 nm Laser 

470 THz with 

-4.6 GHz sidebands 



Microwave ESR 



-4.6 GHz 



Variable Preparation 
Pulse 



Figure Sll. Pulse sequence for arbitrary spin-state readout. The above pulse sequence was used 
for the data presented in Main Text Fig. 4A, B. The DBP readout experiments consisted of 3.75 x 10 6 
iterations of this pulse sequence. 

Our DBP optical spin readout protocol begins with preparation with non-resonant 532 nm ex- 
citation to prepare the state in \0g). This is followed by an on-chip microwave pulse to prepare 
the spin at various points about the Bloch sphere. In the case of Fig. 4A, the pulse corresponds 
to a n/2 pulse, and its phase is varied to place the spin at various points about the Bloch sphere 
equator. In the case of Fig. 4B, its pulse duration is varied to induce Rabi oscillations between the 
ground state spin sublevels, \0g) and |+lg), initializing the spin at various points along a meridian. 
A DBP pulse of 400 ns is used to read out the spin state and is delayed 450 ns from the microwave 
initialization pulse, much like the delay between the CPT initialization and the ESR n/2 projection 
pulses for quantum state tomography of CPT and SRT Spin readout protocols such as QST could 
be performed directly with DBP mitigating the need for ESR; however, the differing PL contrasts 
between polar and azimuthal spin readout as a result of qubit dephasing must be appropriately 
calibrated out (Fig. 4A,B). 

To determine the quality of DBP spin readout vs. traditional "green" spin readout via the ISC, 
we analyze the signal-to-shot-noise ratio of both readout methods. Shot noise of counted photons 
is the primary noise source in our data, which is a common feature for these types of experiments 
with proper mitigation of any systematic errors such as experimental drift. Each data point in 
Fig. 4A,B consists of summed photon counts of n = 3.75 x 10 6 individual spin readouts using a 
400 ns DBP pulse. When fit to a sinusoid, the polar Bloch sphere readout in Fig. 4B results in 
Ibz — 4,850 counts if the spin is in the bright state and 7dz = 1/750 counts if the spin is in the dark 
state (averaged fit values of |0g) and data). The equatorial Bloch sphere readout in Fig. 4B 

with reduced contrast primarily due to dephasing results in I%x — 5,380 counts if the spin is in the 
bright state and iox = 3,160 counts if the spin is in the dark state (averaged fit values of \X g ) and 
\ — Xg) data). To compare to traditional ISC "green" readout, we will assume reasonable numbers 
for green readout with our specific confocal setup: a 400 ns readout window with 20,000 Cts/ s for 
the m s = spin state, and a 30% reduction of counts for m s = ±1 spin states during this window. 
For an equivalent averaging time (n — 3.75 x 10 6 ), this corresponds to green spin readout counts of: 
Ibg — 30,000 and I^q = 21,000. The number of spin readouts (N,) required for the noise standard 
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deviation (y/n) to be roughly equal to the full-scale readout contrast is: 



n /2 x (I Bi + l Di ) 
(hi ~ hi? 



Ni = •li^^l 1 p (S13) 



the number of individual spin readouts to get a signal-to-noise ratio of unity for these three readout 
techniques (Polar DBP (Nz), Azimuthal DBP (Nx), and green ISC (Nq) are respectively: 



N z = 1,290, N x = 3,250, N G = 1,180 

Therefore, using these rough numbers, DBP on a polar axis ("spin up" vs. "spin down") has 
a signal-to-noise ratio comparable to green spin up /down readout, while DBP on a precessing 
equatorial axis ("spin left" vs. "spin right") takes roughly three times as much averaging to get the 
same signal-to-noise ratio as the other two techniques. 
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All-optical Ramsey measurement 



We use the pulse sequence in Fig. S12 to take the data for the top portion of Fig. 4C and in Fig. 
15151 



Traditional 
Initialization 



532 nm Laser 



CPT I nit. & 
DBP Readout 




CPT 




Free precession 

T 




Tunable 637 nm Laser 
- 470 THz with 
DBP -4.6 GHz sidebands 




Figure S12. Pulse sequence for arbitrary spin-state readout. I. We begin with a non-resonant 
532 nm excitation to prepare the state in |0 g ). II. The state is then initialized onto the equator with 
a CPT pulse of length 50 ns, in one of four azimuthal positions n/ 2 out-of-phase from one another 
other. III. The spin state dephases during free precession, x, which is varied. IV. The state is then 
read out with a DBP pulse of length 50 ns of a fixed phase. The data of two opposite CPT pulse 
phases were subtracted and plotted in the top of Fig. 4C. The EOM microwave detuning,(a; mw /2/T — 
Sos/h), is ~7.5 MHz from the mean qubit precession. 



The Ramsey data plotted in Fig. 4C of the main text shows the difference in measured PL for 
two orthogonal initial spin projections, APL = PL( — Xg)) — PL( Xg)). The difference cancels the 
contribution to the PL from the ISC decay in the optical Ramsey measurements (Fig. S13l, and 
effectively measures the projection (Sx) of the spin at the end of the free precession period. In both 
cases the data is fit to a function of the form 

APL = A exp (-^2 J {Cicos[(^o;-a;HF)(T-To)] 



+ cos[Sco(t - t )] + C 2 cos[(Sco + o; hf )(t - r )] j, (S14) 



where 



Scv = w mw - 5 GS /h, (S15) 

and includes the threefold hyperfine coupling to the 14 N nuclear spin with frequency Cl>hf, in- 
homogeneous dephasing with characteristic time TZ, and independent amplitudes for the three 
hyperfine components. The temporal offset Tq accounts for the effects of finite-duration initializa- 
tion and readout pulses. Best-fit parameter values for the curves plotted in Fig. 4C are provided in 
Table 152 
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450 900 1350 1800 

Free precession (ns) 

Figure S13. Non-subtracted all-optical Ramsey measurement. Four different phases of a time- 
domain Ramsey experiment were measured and are plotted above. The overall background is due 
to the CPT pulse populating the dark ISC, which subsequently decays over time, making the DBP 
PL brighter as it gets further away from the CPT pulse. 





i 2 


boo/In 


<^hf/27t 


TO 


A 


Ci 


c 2 




(f/s) 


(MHz) 


(MHz) 


(ns) 


(Cts) 






Optical 


1.13 ±0.05 


7.52 ± 0.01 


2.19 ± 0.01 


13 ±1 


253 ± 13 


1.36 ± 0.07 


0.64 ± 0.05 


ESR 


1.01 ± 0.03 


7.34 ± 0.01 


2.20 ± 0.01 


-12 ±2 


1230 ± 40 


1.16 ± 0.04 


0.68 ± 0.04 



Table S2. Best-fit parameter values from fits of the model of Eq. 
main text. Uncertainties are standard error. 
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to the data in Fig. 4C of the 
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All-optical Hahn echo measurement 



We use the pulse sequence in Fig. S14 to take the data in Fig. 5C. 



Traditional 
Initialization 



CPT Init. & 
DBP Readout 



SRT Rotation 



532 nm Laser 



CPT 



-t/2- 



■t/2- 



Tunable 637 nm Laser 
470 THz with 
DBP -4.6 GHz sidebands 



SRT 




Tunable 637 nm Laser 
470 THz with 
-4.6 GHz sidebands 



Figure S14. Pulse sequence for all-optical Hahn echo measurement. I. A 532 nm excitation 
is used to prepare the state into |0 g ). II. The spin is initialized on the Bloch equator with a CPT 
pulse (Fig. 3A). III. Dephasing of the spin state occurs during free precession for a delay t/2. IV. 
The spin is then rotated by a SRT pulse (Fig. 5A). V. Rephasing of the spin state occurs during a 
period of r/ 2. VI. Finally the spin state is readout along the equator with a DBP pulse (Fig. 4A), 
corresponding to a bright state either in-phase or n out-of-phase with the CPT pulse. The Hahn 
precession time, r, is varied. 



In Fig. 5C, the all-optical Hahn echo measurement is plotted on top of the data set for the ESR- 
based Hahn echo. We use least squares to fit the function APL = Aexp ( — (T/T2) 3 ) to our data 



and infer T2 ~900 \xs for both measurements. The fitting parameters are provided in Table S3 





T 2 


A 




(fiS) 


(Cts) 


Optical 


893 ± 51 


538 ± 29 


ESR 


909 ± 30 


2991 ± 99 



Table S3. Best-fit parameter values from fits of the model to the data in Fig. 5D of the main text. 
Uncertainties are standard error. 
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